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ABSTRACT 



Astrophysical false positives due to stellar eclipsing binaries pose one of the 
greatest challenges to ground-based surveys for transiting Hot Jupiters. We have 
used known properties of multiple star systems and Hot Jupiter systems to pre- 
dict, a priori, the number of such false detections and the number of genuine 
planet detections recovered in two hypothetical but realistic ground-based tran- 
sit surveys targeting fields close to the galactic plane {b ~ 10°): a shallow survey 
covering a magnitude range 10 < V < 13, and a deep survey covering a magni- 
tude range 15 < V < 19. Our results are consistent with the commonly-reported 
experience of false detections outnumbering planet detections by a factor of ~10 
in shallow surveys, while in our synthetic deep survey we find ~l-2 false de- 
tections for every planet detection. We characterize the eclipsing binary con- 
figurations that are most likely to cause false detections and find that they can 
be divided into three main types: (i) two dwarfs undergoing grazing transits, 
(ii) two dwarfs undergoing low-latitude transits in which one component has a 
substantially smaller radius than the other, and (iii) two eclipsing dwarfs blended 
with one or more physically unassociated foreground stars. We also predict that 
a significant fraction of Hot Jupiter detections are blended with the light from 
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other stars, showing that care must be taken to identify the presence of any un- 
resolved neighbors in order to obtain accurate estimates of planetary radii. This 
issue is likely to extend to terrestrial planet candidates in the CoRoT and Kepler 
transit surveys, for which neighbors of much fainter relative brightness will be 
important. 

Subject headings: techniques: photometric — binaries: eclipsing — planetary 
systems 



1. Introduction 



The majority of extrasolar planets that have been discovered by ground-based transit 
surveys to date are so-called "Hot Jupiters", with radii and masses similar to Jupiter and 
orbital periods less than 10 days (d). Ground-based transit surveys can be divided into 
two main types: (i) shallow surveys that use telescopes with ape rtures ^lOc m to monitor 



~50 deg 2 field of views (FO Vs) over a magnitude range V < 13 (eg. lAlonso et al.ll2004 iBakos 
20071 : IStempels et al. I l2007h . and (li) deep su rveys that use ~lm t elescopes to monitor ~1 deg 2 



FOVs over a magnitude range V < 19 (eg. lUdalski et al.ll2003t iBayliss et al.ll2009aU bl) 



Most transit searches have adopted the shallow survey approach, largely because the 
required equipment is less expensive and more readily available than that needed to conduct 
a deep survey. Another appeal of shallow surveys is that ~lm class telescopes can be used 
to perform follow-up spectroscopy on their bright targets, with signal-to-noise ratios (S/N) 
that are sufficiently high to detect radial velocity variations caused by an unseen Jupiter- 
mass companion. The fainter targets of deep surveys, on the other hand, require follow-up 
spectroscopy to be performed on ~6-10m telescopes in order to achieve the same level of 
precision. Certainly, spectroscopic investigations of planetary atmospher es are only currently 
possible for the bright sy stems monitored by the shallow surveys (eg. 
20021 : iKnutson et al.ll2007h . 



Charbonneau et al. 



Nevertheless, shallow surveys do have a number of disadvantages compared to deep 
surveys. Perhaps most significantly, late-type main sequence stars are much rarer over the 
magnitude range covered by shallow surveys than they are in deeper surveys. More stars 
must therefore be monitored overall in a shallow survey than in a deep survey to have the 
same probability of observing a Hot Jupiter transit a late-type main sequence parent star. 

Another drawback of ground-based shallow surveys is the lower spatial resolutions that 
they achieve, typically limited by point-spread functions (PSFs) with full-width-at-half- 
maximum (FWHM) sizes ~20", whereas deep surveys typically have much smaller PSFs 
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of ~2". Spatial resolution is important in transit surveys, because blending with physically 
unassociated stars can dilute any actual planetary signal or cause a transit by an eclipsing 
binary (EcB) system to be mistaken for a planetary system. Although the number of stars 
per square degree of sky is much higher at the faint magnitudes covered by deep surveys, we 
find that the difference in PSF sizes is the dominant effect, making blending more prevalent 
in shallow surveys (Tableland §3.5|) . 

The planetary y ield of any photometric transit survey is far outweighed by the number 
of detected EcBs (eg. IWeldrake fc Baylissll2008t iDevor et al.ll2008l ) . Generally, the light curve 
produced by an EcB is easy to distinguish from one produced by a Hot Jupiter transiting 
its parent star. The transit durations are typically longer and the primary eclipse depths 
are typically deeper for EcBs than they are for transiting Hot Jupiters. Also, because both 
components of an EcB are luminous, a detectable secondary eclipse is usually produced. 
Unless both components of the EcB have very similar radii and luminosities, the secondary 
eclipse will not have the same depth as the primary eclipse, immediately revealing the non- 
planetary nature of the signal. 

While most EcBs can be identified using one or more such diagnostics, in certain cases 
an EcB can mimic the photometric signal produced by a transiting Hot Jupiter. There are 
three main situations in which this can occur: (i) a late-type dwarf star eclipsing a larger, 
brighter dwarf star; (ii) two dwarfs undergoing grazing eclipses; (iii) two dwarfs eclipsing 
each other blended with the light of one or more unresolved stars. Meanwhile, detached 
EcBs wit h a giant component do not exist for periods less than 10 days due to Roche lobe 
overflow (ISetiawan et al.l 12004 ). and hence are not able to mimic the signal of a transiting 
Hot Jupiter. 

It has not been uncommon for the rate of false planetary detections d ue to EcBs to be 
over an order of magnitude higher than the rate of true planet detections (eg. lO'Donovan et al 
20061 ). The identification and rejection of such objects can be far from trivial and poses a 
significant challenge to searches for transiting planets. 

It is instructive to consider the yields of transit surveys in terms of what is already 
known about the properties of mult i ple st ar systems and Hot Jupiter systems from published 
empirical studies. iBeatty fc Gaudil (120081 ) have developed a methodology to predict, a priori, 
the number of planetary detections that a user-specified survey will make. Such estimates 
are helpful both in d esigning surve ys and in evaluating the effectiveness of existing surveys. 
On the other hand, iBrownl ( 120031 ) made a priori estimates of the rate of false detections 
in two synthetic transit surveys, one based on the shallow STARE survey and the other 
based on the then-planned Kepler space-based survey. A bottom-up understanding of the 
types of EcBs that are likely to cause false detections and the rates at which they occur 
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is valuable both for interpreting the objects detected by transit surveys and for refining 
follow-up approaches used to reject non-planetary signals. 



In this p aper, we describe a method analogous to those of iBeatty fc Gaudil ( 120081 ) and 



Brown! ( 120031 ) to estimate the expected rates of false detections and Hot Jupiter detections 
in two synthetic surveys: (i) a ground-based shallow survey; and (ii) a ground-based deep 
survey. For each synthetic survey, the types of EcB configurations that are predicted to 
cause false detections are examined in some detail. These results are then compared with 
the reported experience of real transit surveys surveys, and the general implications for 
ground-based transit surveys are discussed. 

Primarily, our study builds upon that of Brown (2003), by making more rigorous a 
number of the simplifying assumptions contained in that earlier work. For instance, unlike 
Brown, we explore the potential for genuine Hot Jupiter detections to be combined with faint 
blends, which can affect the estimated planetary radius ( §3.31) . whereas he only investigated 
unblended Hot Jupiter detections. We also do not ignore EcB secondary eclipses and use a 
more realistic S/N criterion in our definition of what constitutes a detection ( §3.4p . We have 
also attempted to account for the possibility that blended neighbors fainter than the "formal" 
upper magnitude limit of the survey can contribute in a non-negligible way to the observed 
photometric signals of stars within the formal magnitude range ( §3.5)) . Furthermore, Brown 
only considered single-star blends in his calculations, but we have allowed for multiple-star 



blends ( ^3.6.4p . Our experience with Lupus-TR-3b in the SuperLupus survey (IBayliss et al. 



2009al ibl). which we describe in §21 suggests that such cases might not be especially rare. 



2. Astrophysical false positives 



Survey teams often report that false det ections due to EcBs ou tnumber true planet 
detections by over an order of magnitude (eg. lO 'Donovan et al.l 120061 ). It should be noted 



that reports from shallow surveys dominate the literature, and this ratio may in general be 



lower for deeper surveys, as has been the experience of the SuperLupus survey (IBayliss et al. 



2009aJJb(). In any case, astrophysical false positives are a ubiquitous feature of any transit 
survey and the process of identifying and rejecting them has proven to be a time-consuming 
and, in many cases, subtle task. 

Hot Jupiter transit signals can be imitated by: (i) a late-type main sequence star eclips- 
ing a larger main sequence star, (ii) two main sequence stars undergoing grazing eclipses, 
or (iii) two eclipsing main sequence stars blended with one or more other stars inside a 
photometric confusion radius, determined by either the instrumental PSF or the CCD pixel 



-5 - 



scale. For this reason, survey teams must follow up all promising candidates with additional 
observations to rule out the possibility that the photometric signal is caused by an EcB. 
Such observations include: 



Multi-color photometry: Apart from high-order stellar limb-darkening effects, plane- 
tary transit signals are color independent. Signals due to EcBs, on the other hand, 
typically have a stronger color dependence that can be detected by making photometric 
measurements in more than one passband. 

High spatial resolution imaging: Often a candidate that was originally presumed to be 
a single star will be revealed to consist of two or more fainter stars by higher resolution 
imaging. By monitoring the flux changes of each component separately, it is sometimes 
possible to identify astrophysical false positives caused by blended EcBs by recognizing 
that one of the resolved components undergoes an eclipse that is too deep to be caused 
by a planetary transit. 

Spectral typing: Spectroscopy can be used to identify the spectral type and luminosity 
class of the star that is being transited. If the star turns out to be a giant, then the 
photometric signal cannot be caused by a transiting planet since it would produce a 
flux change too small to be detected. Instead, the signal must be due to a blended 
EcB. 

Radial velocity measurements: For a solar-type star, variations in radial velocity caused 
by an orbiting stellar companion will be ~10kms _1 , whereas the oscillations will only 
be ~0.1kms _1 for a Jupiter-mass companion. Radial velocity measurements with a 
precision of ~1 kins -1 can therefore be used to identify cases in which the companion 
is stellar. 



Latham! ( 120031 ). ICharbonneau et all (120041 ) . and 10 'Donovan et all (120061 . 120071 ) describe 
astrophysical false positives encountered as part of the Vulcan and TrES surveys, detailing 
how the non-planetary natures of candidate objects were revealed using one or more of the 
follow-up techniques listed above. 

An example of a more challenging false positiv e, and the resourc e -inten sive processes 
that were employed to reject it, has been detailed by iMandushev et al.l ( 120051 ). In this case, 
no color dependence was observed in the eclipse and initial spectral typing showed that the 
primary, GSC 01944-02289, was an F dwarf. Radial velocity measurements indicated that 
the companion had a mass of ~30 Mj. All of these observations were consistent with a 
brown dwarf transiting the primary. However, careful measurement of asymmetries in the 
spectral lines and extensive modelling of possible blend scenarios ultimately showed that the 
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observed signal was caused by a faint EcB blended with another unresolved star rather than 
a transiting brown dwarf. 

The experience of the SuperLupus survey faavliss et all l2009al lbl) has also reinforced 



how difficult it can be to determine the precise nature of candidates, especially for the faint 
objects targeted by deep surveys. One promising candidate that has been identified as part 
of this survey, Lupus-TR-3b, appears to be a Hot Jupiter with radius 0.9Rj and mass 0.8Mj 
based on its light curve and follow-up spectroscopy performed at the Magellan II telescope. 
However, high resolution imaging data obtained with the Magellan I telescope, combined 
with image deconvolution analysis, has revealed that the host "star" in fact consists of up 



to se ven separate stars that are unresolved in the SuperLupus photometry (jSackett et al. 



20091 ). While four of these stars are too faint to be responsible for the observed transit 
signal, it is unclear which of the remaining three stars is undergoing eclipse. Therefore, 
while unlikely, Lupus-TR-3b can still not be completely ruled out as being a blended EcB 
rather than a transiting planet. At the very least, if Lupus-TR-3b is a planet, it will have 
different properties to those listed above, which were derived under the assumption that the 
signal was caused by a planet transiting a single star. 

These two examples illustrate the particular challenge posed to transit surveys by blend- 
ing, both in terms of the potential for false detections to be made and the effect that it can 
have on estimated planetary properties. The a posteriori processes used to identify blended 
neighbors, such as spectral line asymmetry analysis and iterative image deconvolution, are 
typically time-consuming and resource-intensive tasks, albeit important ones. Consider- 
ing the problem in the forward direction, starting with well-understood statistics of stellar 
crowding and property distributions of multiple stars and Hot Jupiters provides another 
perspective on the problem, and is the approach taken in this work. 



3. Methodology 

3.1. The synthetic surveys 

We investigate the rates of astrophysical false detections (i.e. EcBs that mimic a Hot 
Jupiter signal) and genuine Hot J upiter detections sep arately for a synthetic deep survey, 
based on the SuperLupus survey (jBavliss et al.l l2009alJbl) . and a synt hetic shallow survey . 



loosely based on surveys such as HAT-Net JBakosI l2007h and TrES jAlonso et al.l booA . 



The important properties of each survey are listed in Table [TJ For simplicity, and ease of 
comparison, the field of views (FOVs) of both synthetic surveys are taken to be centered on 
the same Galactic coordinates as the SuperLupus survey I = 331.5°, b = 11°. The magnitude 
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range of the deep survey is taken to be 15 < V < 19, the same as SuperLupus, while for 
the shallow survey a magnitude range of 10 < V < 13 is used. Using these parameters, 
expected star counts are binned according to luminosity class, spectral type, and apparent V 
magnitude for each synthetic survey using the Besangon model of the Milky Way Galaxy 
(IRobin et al.l I2003T ) . Main sequence stars are found to comprise approximately 75% of all 
stars in our synthetic deep survey, and about 20% in our synthetic shallow survey. 



3.2. Estimating detection rates 

We estimate the rate of false detections Qfd due to EcBs by integrating a probability 
density function Q FD of the form: 

M Q 

$ FD (*,V p ,P,q,cosi,f B ) (1) 



d^> ■ dV p ■ dP ■ dq ■ d(cos i) ■ df B 



where \I/ is the spectral type of the primary, V p is the apparent magnitude of the primary, P is 
the orbital period, q is the secondary-to-primary mass ratio, i is the orbital inclination, and f B 
is the total blended flux. These properties determine the observed primary transit depth 5 P , 
observed secondary transit depth S s , and transit duration r (see §3.6.6p . Hence in effect, 
the integrated density function $_fd gives the probability that any given star monitored will 
be an EcB consisting of two dwarf components with light curve properties {5 P , 5 S , r} that 
satisfy our criteria for a false detection, which are presented in §3.41 A detailed description 
of the integration process is given in §3.61 



The rate of Hot Jupiter detections Sl B j is estimated in exactly the same way, using a 
probability density function of the form: 

d 6 £1 

- $ HJ (*,V p ,P,R,coai,f B ), (2) 



■ dV p ■ dP ■ dR ■ d(cos i) ■ df B 
where the binary mass ratio q is replaced with the planetary radius R 



3.3. EcB and Planet-Star configurations considered 



Ground-based transit surveys are overwhelmingly sensitive to transiting planets with 
periods P < 10 days, due to observational window effects as well as the low geometric prob- 
ability of a transit occurring for longer periods. Therefore, we only consider EcBs and Hot 
J up iters with periods in this range. This allows us to ignore EcB s with a giant compo nent, 
as detached giant binaries do not exist for periods P < 10 days (ISetiawan et al.l 12004 ). We 
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also only consider Hot Jupiters transiting main sequence stars, since the p hotometric signal 



cause d by a Hot Jupiter transiting a giant star is too faint to be detected (IGould fc Morgan 



20031 ). 



We divide both false detections and Hot Jupiter detections into two main types: blended 
and unblended. Blended detections are those for which the eclipsing system is accompanied 
by one or more unresolved neighbor stars. The resolution limit r c is taken to be approximately 
equal to the size of the PSF or the pixel scale, whichever is largest. For the deep survey 
we use a value of 1.5", which is roughly equal to the spatial resolution achieved by the 
SuperLupus survey, and for the shal low survey we use a value of 20", which is typical of the 



PSFs in actual shallow surveys (eg. iBakos et al.ll2004l ). Blending has the effect of reducing 



the measured transit depths 5 m , according to: 



5 m = f {" f K (3) 

Ju + JB 

where /„ is the flux of the unblended EcB, f B is the total flux of the blended neighbors, and 
S u is the transit depth that would be measured if the EcB was unblended. 

With regard to tertiary companions of EcBs, in reality some would fall within the con- 
fusion radius r c and some would fall outside. However, we chose not to estimate the true 
fractions of resolved and unresolved tertiary components, and instead perform the calcula- 
tions separately under two limiting assumptions: (i) all tertiary companions are resolved, 
i.e. fall outside the confusion radius r c of the survey photometry; and (ii) all tertiary com- 
panions are unresolved. In the latter case, the tertiary companion is treated as an additional 
blended flux in Eq. [3J Higher order companions are not accounted for. 

We do not take stellar multiplicity into account when calculating Hot Jupiter detection 
rates. That is, we treat every star as single and ignore possible complications associated with 
detecting Hot Jupiter transits in multiple systems such as blending by binary companions. 



Beatty fc Gaudil (120081 ) also made this simplification, and estimated that it could result in an 
overestimation of planet detection rates by at most 30%. However, we do consider blended 
cases of Hot Jupiter detections (i.e. planet + host star + physically unassociated stars), in 
addition to unblended Hot Jupiter detections. 



3.4. Criteria for a detection 



We define an astrophysical false positive as an EcB that meets all of the following 
criteria: 
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Both stellar components are dwarfs with masses greater th an O.IMm, which approxi - 
mately corresponds to the minimum hydrogen burning mass Jchabrier fc Baraffelboool ). 



We do not include transits by brown dwarf companions in our false detection counts, 
since brown dwarfs orbiting within 5AU of main sequence star s are known to be ex- 



tremely rare (the "brown dwarf desert") (IMarcy fc Butlerl 120001 ) . 



The observed primary transit depth is 5 P < 0.05 and has S/N> 7, where 7 is some 
minimum detection threshold. While it is possible for a Hot Jupiter transiting a late 
M dwarf to cause a fractional flux change greater than 0.05, such cases are rare, and 
generally disregarded by survey teams due to the high probability that they are caused 
by an EcB. 

The observed secondary transit depth 8 S either has S/N< 7 or 5 P — S s < e, where e is the 
measurement noise, i.e. it is either too shallow to be detected or it is indistinguishable 
from the primary eclipse. 

The transit duration r is less than 0.25 d, which approximately corresponds to the time 
taken for a Hot Jupiter with a period less than 10 d to undergo a central transit of the 
disk of an F dwarf. Hot Jupiter transits of dwarfs earlier than F, which have larger 
radii and hence longer transit durations, will g enerally be too shallow for detection by 



ground-based surveys (IGould fc Morgan! 120031 ) . EcBs with periods P < 10 d, on the 



other hand, are much more likely to have transit durations r > 0.25 d due to the large 
relative disk sizes of the eclipsing bodies. 

5. The effective magnitude V e ff of the EcB combined with any blended neighbors is within 
the magnitude limits of the survey. 



With the exception of the mass constraint on the secondary, our criteria for a Hot 
Jupiter detection are exactly the same as those listed above for an EcB false detection. Of 
course, the secondary eclipse that occurs when the planet passes behind the star does not 
produce an observable signal, so for a transiting Hot Jupiter condition (3) will always be 
satisfied. We also note that Hot Jupiter transits with primary transit depths 5 P > 0.05 or 
transit durations r > 0.25 d are possible, but since such signals are far more likely to be 
caused by EcBs, survey teams generally discard such candidates if they arise. 



3.5. The extended magnitude range 



In addition to the "formal" survey magnitude range, which we denote by V& < V < Vf 
(i.e. 10 < V < 13 for our synthetic shallow survey, and 15 < V < 19 for our synthetic deep 
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survey), we define an extended magnitude range V, < V < V e for each survey, where Vf <V e . 
This is to allow for the possibility that blended stars with apparent magnitudes fainter than 
the formal survey cut-off Vf can contribute to the observed photometric signal. The faint 
limit V e of the extended magnitude range is taken to correspond to the flux f e ~ a m , where 
a m is the lowest root-mean square (rms) scatter a in the photometry achieved over the formal 
magnitude range. That is, as long as a star has flux at least at the level of the noise, we 
account for its contribution in the total measured flux. 

To obtain analytic expressions for the fractional rms a as a function of apparent mag- 
nitude V , we fit parabolic curves of the form: 

log 10 a = aV 2 + bV + c (4) 

by eye to the scatte r in light curves as a function of V in the Supe rLupus survey (Fig. 2, 



Bayliss et al.ll2009bl ) and the TrES survey (Fig. 2.2. IO'Donovanll2007l ). corresponding to our 
synthetic deep and shallow surveys, respectively. Values for a, b, and c for each survey are 
listed in Table [TJ Using Eq. HI we find V e = 22 for the deep survey and V e = 16 for the 
shallow survey. 

In order for any star with magnitude Vf < V < V e to contribute to a photometric 
measurement, however, it must be blended with at least one other star that has magnitude 
V, < V < Vf, i.e. within the formal survey magnitude limits. To simplify the analysis, 
we ignore those cases in which two or more unresolved stars each have a brightness fainter 
than the faint limit of the survey, but a total brightness greater than the faint limit. Such 
"borderline" cases will be relatively rare and should not have a significant effect on the final 
calculated result. 

Supposing we have a FOV containing a total of N stars, the probability p n that any one 
of those stars will be blended with precisely n other stars is given by 

(AT _ 1)1 / ^JV-n-l 



Vn n\(N-n-l)\M" [} Ad) (5) 

where M = Apov/^^ is the number of resolution elements in a FOV of area Apovy given 
a confusion radius r c . 

By letting N in Eq. [5] be equal to the total number of stars a in the FOV (both resolved 
and unresolved) with magnitude < V < Vf, we obtain the fraction p n of stars that are 
blended with exactly n other stars with magnitudes in the formal magnitude range. This 
means that if /3 is the number of stars with magnitude Vf < V < V e between the faint limit 
Vf of the formal magnitude range and the faint limit V e of the extended magnitude range, 
then only (1 —po)/3 of them will be able to contribute to the photometric measurements, i.e. 
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those that are blended with at least one other star within magnitude in the formal range. 
The effective total star count over the extended magnitude range is therefore: 

iV s = a+(l-p )/3 (6) 

and the number of resolved objects Nr is given by: 

n 

In §4.1[ we use Eq. [7] to quote our detection rates as the number of detections per N R = 
10, 000 resolved stars. 

Values for the crowding probability p n in each synthetic survey computed using Eq. [5] 
and |6] are listed in Table EJ The low blending probabilities p n for n > 1 we obtain for the 
synthetic deep survey could mean that the case of Lupus-TR-3b, which consists of up to 
seven blended stars (§|2]), is an unlikely single event. On the other hand, of the seven blended 
stars, at least four have magnitudes V > 22, making them too faint to be accounted for in 
our analysis which imposes a str ict cut-off at V = 2 2. Two of the remaining blends may also 



be a wide physical binary pair ( iSackett et al.l 120091 ) . It is also possible that the star counts 



we obtain from the Besangon model are underestimated. For instance, the actual interstellar 
extinction may be lower than the value of O.Tmagkpc^ 1 that we use for the Besangon model 
input. 



3.6. Integrating the probability density function 

To estimate detection rates, we integrate the probability density functions given by Eq. [1] 
and [2] over appropriate property ranges (see below). Due to the independence of parameters 
in these relations, the probability density functions &fd and can be broken into a 
number of component probabilities. For the false detections: 

$ FD (V,V P , P, q, cos i,f B ) = <S> d (y,V p )-$ b (P,q)-$i(cosi)-$ f (f B ) 

x $ W (P) .<$> lc {V,V p ,P,q, cos i,f B ) (8) 

where V p ) is the probability that the star will be main sequence with spectral type \l/ 

and apparent magnitude V p ; q) is the probability that a main sequence star will be 

the primary component of a multiple system with period P and mass ratio q; $j(cosi) is 
the probability that the orbital inclination is %\ is the probability that the total 

flux due to blended companions is fs] <& W (P) is the probability that at least three tran- 
sits will be observed given the observing window function; and $/ c (\l/, V p , P, q, cosz, is 
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the probability that the transit signal of a binary consisting of two dwarfs with proper- 
ties {\&, V p , P, q, cosi, f B } resembles that of a transiting Hot Jupiter. The assumption of 
independence between the spectral type \l/ and the other propert ies of the binary, nam ely, 
the period P and mass ratio q, is supported by empirical studies (IHalbwachs et al.ll2003l ). 



Similarly, for the planet detections, we can divide the probability density function of 
Eq. 12 into: 



$ H j(V,V p ,P,R,cosi,f B ) 



M*,V P ) ■ $ fe (P,i?) ■ ^(cosz) • <f> f (f B ) 
x $ w (P)-<f> lc (V,V p ,P,R,cosi,f B ) 



(9) 



where $&(P, R) refers to the probability that a main sequence star will host a Hot Jupiter 
companion with orbital period P and radius R, and all other terms are the same as for the 
false detection case given by Eq. [SJ Again, we assume independence between the spectral 
type of the primary and the companion properties, which in this case are the planet 
radius R and orbital period P. While correlations are beginning to emerge between system 
parameters as m ore transiting planet systems are discovered, such as between stellar m ass 
and planet mass (IJohnson et al.ll2007l ). and planet mass and period (ITorres et al.ll2008l ). we 
note that the observed scatter is still larg^]. Therefore, while our assumption of independence 
between the system parameters mentioned above is a simplification, it is reasonable for the 
purpose of this study and is unlikely to have a substantial effect on our final results. 

For the EcB false detections, integration is performed numerically over the following 
ranges: 



• All spectral types ^ G [OO, M9], with a step size d^f of five spectral sub-types. 

• The extended magnitude range V p G [15, 22] for the deep survey and V p G [10, 16] for 
the shallow survey, with a step size dV p = 0.5 mag. 

• Periods P G [1, 10], with a step size dP = 0.5 days. 

• Mass ratios g 6 [0, 1] with a step size dq = 0.1. 

• Orbital inclinations cosz G [0, 1] with a step size d(cosi) = 0.01. 

• For the blended cases, fluxes f B G [//, /&] where ff is the flux corresponding to the 
faint magnitude limit of the extended magnitude range and fb is the flux corresponding 
to the bright magnitude limit. A step size df B corresponding to half a magnitude is 
used. For the unblended cases, f B is not integrated over, i.e. $/(0) = 1 in Eq. [8]and[9j 
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For the Hot Jupiter detections, spectral types \& G [FO, M9] are integrated over, since 
dwarfs of s pectral type earlier than about FO are too large for a Hot Jupiter transit to be 
detectable ( IGould k, Morgan! 120031 ). We integrate over Hot Jupiter radii R G [0.8, 1.6]Rj 
with a step size dR = O.lRj. We choose this range of radii because it encompasses the vast 
majority of Hot Jupiters discovered via transits at the time of writing^. All other integration 
ranges are the same as for the EcB false detections listed above. 

In the following subsections, we describe the observational data and assumptions that 
are used to define the individual components of the probability densities in Eq. [S] and EO 



3.6.1. Stellar property distribution ^ d (^,V p ) 



We use star counts from the Besangon model (IRobin et al.ll2003l ) to calculate the proba- 
bility V p ) d^> dV p that a star will be main sequence with spectral type \1/ and apparent 
magnitude V p : 



N(*,V P . 



(10) 



where N(^>, V p ) is the number of main sequence stars with spectral type \l/ and apparent 
magnitude V p , and iV 2 is the number of stars of all luminosity classes, spectral types, and 
magnitudes in the extended magnitude range given by Eq. El 



3.6.2. Companion property distributions q) and &t,(P,R) 

We assume that the period P of a stellar binary is independent of the mass ratio q, i.e. 



$ b (P,q) = $ P (P) • %(q) 



(11) 



where <&p(P) is the probability that a main sequence star will be the primary of a binary 
system of period P, and <& q (q) is the probability that a main sequence binary system will 
have mass ratio q. The distr ibutions for <&(P) and & q (q) are taken from the empirical study 
of spectroscopic binaries by lHalbwachs et al.l ( 120031 ) . who found no evidence for the period 
being dependent on the mass ratio over the range < P < 10 d. This is consistent with 
our assumption of independence between these two properties (Eq. [TT]) . When performing 
calculations for the limiting case of all tertiary companions being unresolved (see §3.3p . we 
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take the fraction of binary systems with tertiary components as a function binary period 
from iTokovinin et al.l ( 120061 ) . 



Similarly, we assume that the period P and radius R of a Hot Jupiter are independent, 



i.e. 



(12) 



where $p(P) is the probability that an F0-M9 main sequence star will host a Hot Jupiter 
of period P, and $r(R) is the probability that a Hot Jupiter wil l have radius R. Th e dis- 
tribution for $p(P) is obtained from the empirical estimates of iFressin et al.l ( 120071 ). and 
for simplicity we assume a uniform distribution for $p(i?) in the range of radi i we consider 
here, namely 0.8-1.6 Rj. We chose to take the period distribution $p(P) of IFressin et al. 



(120071 ). which was derived from the results of the OGLE tran sit survey, as opposed to anal- 



ogous distributions obtained from radial velocity surve ys (eg. ICumming et al.l 20081 ) due to 
the inherent metallicity bias in the latter. See §2.2 of iBeatty & Gaudil (120081 ) for a fuller 
discussion of this issue. 



3.6.3. Inclination distribution $j(cosz) 
Orbits are assumed to be orientated randomly in three-dimensional space, giving $j(cos i) = 

1. 



3.6.4- Blending distribution $/(/b) 

The probability $/(/b) that a star is blended with one or more other stars with a total 
blended flux f B is 

*/(/*) = (13) 

n 

where 3>/(n, /b) is the probability that the blend consists of exactly n other stars with total 
flux f B , and is given by 

where if denotes all combinations of fluxes fi, f'2, ■ ■ ■ , f n that satisfy 

n 

fB = J> (15) 



i=l 



M n 



1 - 



1 \ 



M 



JV s -n-l 



(14) 
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Ni is the number of stars in the field with flux /j, M is the same as in Eq. [5], and N-£ is the 
effective total number of stars in the FOV as given by Eq. [6j Values for the blended flux 
probabilities computed using Eq. [T3] and [0] are plotted for each synthetic survey in 

Fig.m 



3. 6. 5. Window function Q w (P) 



For both surveys, we use the window function <$> W (P) provided in Fig. 2 of iBayliss et al. 



( 12009al ). This function is an approximation of the real window function for the SuperLupus 
survey, and is calculated assuming 100 contiguous observing nights with weather statistics 
that match actual weather logs from the Siding Spring Observatory site. 



3.6.6. Light curve detection probabilities $/ c (^, V p , P, q, cosi, fs) and 

®ic(^, V p , P, R, cosi, / B ) 

Assuming that the orbit is circular, and given the properties V p , P, q, cosi, we 
can determine the observed transit depths, S p and 5 S , and the transit duration r as follows. 
First, we use the spectral type \I/ of the primary to find its mass M p , rad ius R P: and lumi- 
nosity L p using a fit to the relations for main sequence stars tabulated by I Cox! ( 120001 ) . The 
luminosity combined with the apparent magnitude V p then provides the distance d to the 
system. The primary mass and the mass ratio q yield the secondary mass M s , and hence the 
secondary radius R s and luminosity L s . The distance d and secondary luminosity L s is next 
used to calculate the flux f s from the secondary. The masses, radii, and period P provide 
the orbital semi-major axis a via Kepler's Third Law. 

With these properties, we have enough information to compute the observed primary 
transit de pths, S p and <L , and the transit duration r using standard relations such as those 
set out in [Sackett (1999|). We assume a quadratic limb darkening law with solar coefficients 
taken from Cox ~(j200ol ) . regardless of spectral type. 



The S/N of the transit signature in each light curve is then calculated for each of S p and 



5 S via: 



S/N = ^ (l(i) 
a 

where a is the measured out-of-transit photometric rms scatter given by Eq. HI and \i is the 
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number of in-transit data points, which we take to be: 

rjhv r 



(17) 



where 77 is the number of observing nights, h is the observing time per night, v is the fraction 
of usable nights, and u is the measurement cadence. For simplicity, and ease of comparison, 
we use i] = 100 nights, h — 8 hours, v = 0.67, and uj = 6 minutes for both synt hetic surveys, 
which is comparable with the actual experience of the SuperLupus survey (jBayliss et al. 
2009al lbh. 



Equipped with the transit depths 5 P and 5 S , the associated S/N ratios, the transit 
measurement noise e = a / -JJL and the transit duration, we are then able to determine if the 
EcB configuration satisfies the criteria for a detection given in §3.41 The probability is then 
simply: 



<$>i c {^,V p , P, q, cos i,f B ) 



1 if the detection criteria are met 



otherwise 



The process for the Hot Jupiter detections is directly analogous to that outlined above for 
false detections. The only difference is that instead of a mass ratio q being used to obtain the 
secondary mass, the Hot Jupiter mass M is obtained from the radius R by assuming all Hot 
Jupiters have the same density pj as Jupiter, i.e. M = 47ri? 3 pj/3. This is a simplification, 
considering the observed range in mean densities of the transiting planets discovered so far, 
from abno rmally-inflated planets such as TrES-4 with only ^0.1 times the mean de nsity 
of Jupiter (IMandushev et al.l 1200 7l) . to objects such as XO-3b (IJohns-Krull et al.l 120081 ) and 
WASP-18b f Hellier et al.l 2009 ) with mean densities ~10 times that of Jupiter. However, 
this simplifying assumption will only have a minor effect on the final results. This is because 
we only use the mass of the planet in the calculation of the orbital semi-major axis a, which 
is then used with the orbital inclination cosz to calculate the transit latitude a cost. Since 

1 /3 

the semi-major axis is a ~ (M p + M s ) 1/0 for a given orbital period P according to Kepler's 
Third Law, and M s << M p in the case of a star-planet system, the uncertainty in the mass 
of the planet will only have a negligible effect on the calculated semi-major axis and hence 
on the calculated transit latitude. 

It must be noted that Eq. [TS] assumes Poisson statistics for the noise, i.e. that the noise 
decreases according to ~ 1 /%//•* (white noise). In practice, however, effects such as varying 
airmass, weather conditions, moon phase, and telescope tracking cause systematic trends in 
the measure d light curves th at can be correlated over the timescale of a planetary transit 
(red noise) (IPont et al.l 120061 ) . This means that a measured light curve will in fact have 
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a lower S/N than predicted by Eq. [161 Rather than at tempting to model the white noise 
and red noise components individually, as was done by iBeatty &: Gaudil (120081 ). we simply 
require that the underestimated S/N calculated using Eq. [TH] exceeds a higher cut-off than 
a survey team would typically require. Specifically, we calculate separate detection rates for 
S/N cut-offs of 10, 20, and 30. By varying the S/N cut-off, we are also able to investigate 
the sensitivity of the yield to this component of the detection criteria. 



4. Results 

4.1. Predicted rates 

Calculated detection rates for the synthetic deep and shallow surveys obtained using S /N 
detection cut-offs of 10, 20, and 30 are reported in Tables [3] and HI assuming separately that 
all tertiary components in multiple star systems are resolved and unresolved, respectively 
( §3.3p . Values are expressed as the number of detections per 10, 000 resolved stars monitored 
with sufficient photometric precision to detect a Hot Jupiter transit. Recall that the number 
of resolved stars N R in the survey field is given by Eq. [TJ Values without parentheses 
are those that have been calculated with a window function ( §3.6.5p . while those values 
with parentheses have been calculated ignoring window effects. In both synthetic surveys, 
the inclusion of our window function does not significantly reduce the predicted number of 
detections. This reflects the fact that the window function we used ( §3.6.5p gives almost 
complete recovery for signals with periods < 5 days. 

The results in Tables |3] and H] show that the total false detection rates are hardly changed 
if we assume that all tertiary components are resolved or unresolved in the survey photometry. 
This is not surprising, considering that in our calculations, an unresolved tertiary component 
only has the effect of increasing the overall flux by an amount less than or equal to that of 
the primary (since we assume a mass ratio < 1, §3.3p . Therefore, some of the EcB 
transits that are too deep to be mistaken for a planetary signal if the diluting effect of 
the tertiary component is not accounted for will become shallow enough to cause a false 
detection when the additional flux is included (Eq. [3]). On the other hand, some of the EcB 
transits that are only just deep enough to be detected when they are not blended with an 
unresolved tertiary component will instead fall below the noise when the additional flux of an 
unresolved tertiary component is included in the transit depth calculation. The fact that our 
overall false detection rates show hardly any change if we assume all tertiary components are 
resolved or unresolved indicates that these competing effects are roughly equal in magnitude. 

The results in Table H] suggest that hierarchical triple systems are likely to constitute a 
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substantial fraction of false detections in transit surveys close to the galactic plane. For the 
limiting case that all tertiary components are unresolved, in both synthetic surveys we find 
that hierarchical triples cause ~5 times more false detections than EcBs without tertiary 
components. This result can be exp lained by the fact tha t short period binaries are highly 



likely to have tertiary components ( jTokovinin et all 120061). and the y also have a relatively 



high geometric probability of transiting (see Eq. 2 of ISackettl I19991 ) . 



Overall, we find ~2-4 false detections and ~l-3 Hot Jupiter detections per 10, 000 stars 
monitored in the deep synthetic survey, depending on the S/N cut-off used to define a detec- 
tion ( §3.6.6|) . False detections are divided roughly evenly between blended and unblended 
detections, whereas only about 20% of planet detections are blended. 

These results compare well with the SuperLupus survey in which a total of 13 planetary 
candidates have been identified from a sample of ~25,000 stars monitored (Bayliss 2009, 
private communication). Of these, however, at least five do not meet the all of the detection 
criteria used in the present study. They either have transit durations longer than 0.25 days 
or secondary eclipses that were detected once greater phase coverage was obtained for the 
orbit. This leaves a total of eight candidates, six of which have not yet been followed up. Of 
the two candidates that have been followed up, one is most likely a transiting Hot Jupiter 
(Lupus-TR-3b, $2]) and the other has been confirmed as a false positive. 

In the synthetic shallow survey, we find ~3-4 false detections, similar to the synthetic 
deep survey, but only ~0. 1-0.6 Hot Jupiter detections per 10,000 stars monitored. In other 
words, one planet is detected for every ~15, 000-100, 000 stars monitored, and for every 
planet detected we obtain ~5-25 false detections, depending on the S/N cut-off used in the 
calculations. Approximately 90% of the false detections in the shallow survey are caused 
by blended EcBs, whereas blended and unblended planet detections occur in roughly equal 
amounts. 

These results seem broadly consistent with the commonly-reported experience of shallow 
surveys that astrophysical false positives outnum ber planet detectio ns by an order of mag- 
nitude. For instance, as part of the TrES survey, lO'Donovan J2007h used the 10cm Sleuth 



telescope to monitor ~150,000 stars across 19 fields with photometric precision < 2%. From 
this, 67 candidates were identified but only 4 were confirmed as genuine planet detections. 
By comparison, our results predict for every 150, 000 stars monitored there will be ~2-12 
Hot Jupiter detections and ~40-70 false detections, depending on the S/N cut-off used and 
whether or not a window function is included. The agreement to within a factor of ~2 is 
reasonable, considering the idealized nature of our detection criteria, our assumption that 
full orbital phase coverage is obtained, and the fact that our calculations were made for a 
single field location whereas Sleuth monitored several fields spread over a range of galactic 
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latitudes. 

Another point worth noting about the results for the synthetic shallow survey listed 
in Tables [3] and @] is that the rate of false detections caused by unblended configurations 
remains constant within the numerical precision as the S/N cut-off is increased from 10 to 
30. This is due to the fact that as the S/N cut-off is increased, less primary transits and less 
secondary transits are detected above the noise. The former has the effect of decreasing the 
false detection rate, whereas the latter has the effect of increasing it. Since no net change is 
observed in the calculated detection rates, both opposing effects must be equal in magnitude 
for the unblended false detections in the shallow survey. 

More specifically, the type of configuration that is likely to cause a false detection at 
low S/N cut-offs but not at higher S/N cut-offs are those for which: (i) the primary transit 
signal is only just deep enough to be detected above the noise at the lower S/N cut-offs but 
not at the higher cut-offs; and (ii) the secondary transit signal is either shallow enough to 
escape detection or is indistinguishable from the primary transit. Such configurations will 
lead to a reduction in false detections as the S/N cut-off is increased. At the same though, 
there will be configurations for which the primary transit is deep enough to be readily 
detectable (i.e. S/N> 30) but not so deep that it is ruled out as a planetary candidate (i.e. 
5 P < 0.05), while the secondary transit is only shallow enough to evade detection at the 
higher SN cut-offs. Such cases will increasingly contribute to the number of false detections 
as the S/N cut-off is increased. Both of the situations described here can arise from grazing 
configurations of an F or G star and a smaller star such as a late K or M dwarf (see §4.2p . 
the difference being that the transit latitudes will be lower in the former case than they are 
in the latter. 

In the synthetic deep survey, however, the unblended false detection rate decreases as 
the S/N cut-off is increased from 10 to 30. The reason for this is the lower photometric 
precision attained in the deep survey compared to the shallow survey (Eq. @]and Tabled]). 
As a result, there are relatively fewer configurations with a distinct secondary transit that 
is detectable at lower S/N cut-offs but not detectable at higher S/N cut-offs than in the 
shallow survey. Instead, most of the unblended false detections in the deep survey have 
primary transit depths that are only just detectable above the noise at the lower S/N cut- 
off, and secondary transits that are either too shallow to be detected or of a similar depth 
to the primary transit, as in the first case described above for the synthetic shallow survey. 

For both synthetic surveys, the rate of blended false detections decreases as the S/N 
cut-off is increased. Again, this is due to configurations with primary transits that are 
detectable at the low S/N cut-offs becoming undetectable as the cut-off is increased, and a 
secondary transit that is either too shallow to be detected above the noise or not discernible 
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from the primary transit. This mainly occurs for cases in which significantly more flux comes 
from the blended stars than from the EcB itself ( §4.21) . resulting in shallower transit depths. 
Meanwhile, since planetary occultations only exhibit a primary transit signal, the rates of 
blended and unblended planet detections also decrease as the S/N cut-off is increased in 
both surveys. 

We emphasize that our current study is intended to be illustrative: it should not be 
assumed that the results reported here for a FOV centered on galactic latitude b = 11° 
also apply to surveys conducted further away from the galactic plane, where the stellar 
population statistics and crowding probabilities will be different. A proper investigation of 
how detection rates are affected by the galactic latitude of the survey field is left to a future 
analysis. 

4.2. Predicted detection trends 

The rate of false detections and Hot Jupiter detections in both the synthetic shallow and 
deep surveys are shown in Fig. [2] as a function of the spectral type of the primary star being 
eclipsed. These distributions closely follow the star counts returned by the Besangon model. 
However, in both surveys, the planet detections are skewed towards late-type dwarfs. This is 
due to the deeper, and hence more readily detectable, eclipses that result when Hot Jupiters 
transit late-type dwarfs which have comparatively small radii and low luminosities. For 
the false detections, we found that increasing the S/N cut-off between 10 and 30 biases the 
relative fraction of detections slightly towards earlier spectral types (F5 and GO). This occurs 
because the photometric rms scatter decreases towards brighter magnitudes ( §3. 5j) where the 
fraction of earlier spectral types increases relative to later spectral types. However, this effect 
is quite small, with the relative fractions of detections changing by < 5% for a given spectral 
type, and so we have only included the plotted results for the S/N> 20 cut-off scenario in 

Fig. m 

Fig. [3] and H] show the rate of false detections and Hot Jupiter detections as a function of 
effective magnitude in the deep synthetic survey and shallow synthetic survey, respectively. 
Results are shown for the S/N> 10 and S/N> 30 cases to illustrate the effect of that varying 
the S/N cut-off has on the relative fractions of detections. In both synthetic surveys for the 
S/N> 10 cut-off and S/N> 30 cut-off, the rate of false detections and Hot Jupiter detections 
initially increases with increasing magnitude, simply reflecting the increasing numbers of 
stars at fainter magnitudes. A maximum is then reached, followed by a decline in detection 
rates over the faintest portion of the survey magnitude range. The decline in rates occurs due 
to the decrease in photometric precision at fainter magnitudes. For the same reason, Fig. [3] 
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and H] show that the peak in the relative fraction of detections is shifted towards successively 
brighter magnitudes as the S/N cut-off used to define a detection is increased from 10 to 
30. An interesting feature of Fig. [3] and S] is that the peak in the relative fraction of planet 
detection rates tends to occur at brighter magnitudes with higher photometric precisions 
than it does for false detections, especially as the S/N detection cut-off is increased. This 
is because the number of EcBs with S p ~ 0.05 is greater than for Hot Jupiter transits, with 
the latter tending to produce shallower eclipses that are harder to detect above the noise. 

Fig. [5] shows how the number of blended detections varies with the fraction of blended 
light in both the synthetic deep and shallow surveys. We found that increasing the S/N 
detection cut-off slightly shifted the relative fraction of false detections towards fainter blends 
(i.e. blends contributing less of the overall flux) and the relative fraction of Hot Jupiter 
detections towards brighter blends. However, as in the case of relative detection rates versus 
primary spectral type (Fig. |2J), this effect is small (< 5%), so we have only plotted the results 
for S/N> 20 in Fig. |5j These results show that around 35% of the blended false detections 
in the deep survey are due to bright blends; namely, those blended false detections in which 
the EcB contributes < 30% of the total flux. In the shallow survey, roughly 80% are due to 
bright blends. For blended planet detections in both surveys we find that blends of up to 
about 50% can occur, which has implications for planetary parameter estimates (see §5|). 

The EcB configurations responsible for unblended and blended false detections are il- 
lustrated in Fig. Eland respectively. The primary purpose of these plots is to qualitatively 
illustrate the most common false detection configurations. Results are only shown for the 
case of S/N> 20 since no appreciable change is observed as the S/N detection cut-off is 
increased from 10 to 30. We see that the majority of unblended false detections are caused 
by grazing "twins", which have indistinguishable primary and secondary eclipses. The re- 
mainder of cases are made up of: (i) grazing transits by stars with different radii; and (ii) 
low-latitude transits of a primary by a much smaller companion. For the blended false de- 
tections, there are two broad categories: (i) bright blends with a low- latitude EcB; and (ii) 
fainter blends with a grazing EcB. 

5. Discussion and Implications 

The estimated false detection rates listed in Tables [3] and @] are similar for the deep and 
shallow synthetic surveys, but the estimated planet detection rates are ~5 times higher for 
the deep survey than they are in the shallow survey. There are two main reasons for this: 
(i) ~75% of all stars in the deep survey FOV are dwarfs, whereas only ~20% are dwarfs in 
the shallow survey FOV, according to the Besangon model star counts; and (ii) it is more 
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difficult to detect a blended transiting Hot Jupiter than an unblended one, and the blend 
probability is higher in the shallow survey due to the larger confusion radius (Table [2]) . In 
practice, however, the comparatively low ratio of false detections to planet detections in 
the deep survey is at least partially offset by the inevitable difficulty of conducting suitable 
follow-up observations to verify or reject their faint candidates. 

Extrapolating from our synthetic surveys, ~20% of all false detections in shallow surveys 
and ~65% of all false detections in deep surveys close to the galactic plane will be due to 
EcBs that are either unblended or which contribute more than 30% of the light in a blended 
configuration. For the bright targets of shallow surveys, such cases will generally display 
radial velocity variations that are readily detectable with lm-class telescopes. For the faint 
targets of deep surveys, however, acquiring follow-up spectroscopy with sufficient S/N will 
always be a challenge, but these configurations present the most accessible configurations to 
rule out. 

We find that ~65% of all false detections in shallow surveys close to the galactic plane 
will be caused by bright blends with a physically unassociated giant. The equivalent fraction 
is ~10% in deep surveys close to the galactic plane, reflecting the lower frequency of giants 
over fainter magnitude ranges. Typically, such cases are relatively easy to rule out, either 
through spectral typing or reference to existing stellar catalogs, at least in the case of the 
bright shallow survey targets. 

The other bright blend false detections, which constitute ~15% of all false detections in 
shallow surveys and ~ 25% of all false detections in deep surveys close to the galactic plane, 
will be due to blending with physically unassociated dwarfs. It is this latter case of blended 
EcBs that pose the greatest challenge to follow-up observations. Multi-color photometry 
can often be useful to reveal a color-dependence in the eclipse depth. However, the task of 
unmasking such false positives would be considerably more difficult if the colors of the EcB 
components happen to be similar and the blended dwarf is a rapid rotator, in which case its 
broadened spectral lines may conceal those of the fainter EcB. 

6. Conclusion 

We have investigated the rates of both false detections caused by EcBs and genuine 
Hot Jupiter detections in ground-based transit surveys from a bottom-up perspective, us- 
ing empirically-determined property distributions of multiple star systems and Hot Jupiter 
systems as our input. In particular, we find that in both deep and shallow surveys: (i) signif- 
icant numbers of false detections can be caused by faint EcBs blended with the light of one 
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or more foreground dwarfs, a configuration which experience shows can produce particularly 
insidious impostors, and (ii) for genuine Hot Jupiter detections, up to ~50% of the system's 
light can be due to unassociated blends, which will affect the estimated planetary parameter s 
if they are not accounted for. We therefore reiterate the conclusion of ISackett et al.l ( 120091 ): 
namely, that it is important to identify or rule out blended neighbors for all candidates, 
including blended neighbors that are up to ~5 magnitudes fainter than the primary. 



Lastly, we note that the space-based missions CoRoT (jBaglin et al.l 120071 ) and Kepler 
( jBorucki et al.ll2003l ) may be particularly susceptible to false detections or planet parameter 
misestimation caused by blending. Due to the high photometric precision of these instru- 
ments, blended neighbors up to 8 to 10 magnitudes fainter than the primary can make 
non-negligible contributions to the measured light curves if the signal that is being sought is 
that of an Earth-sized object transiting a dwarf star. Furthermore, since these surveys target 
planets with masses down to about an Earth mass, and the monitored stars are quite faint 
( ^ ~ 15-16), radial velocity follow-up will be of little use in confirming their low-mass can- 
didates. Constraining the presence of faint blended neighbors will therefore be of particular 
importance if a terrestrial-sized planet is to be secured with confidence in these surveys. 



T.M.E. would like to thank The University of Sydney Department of Physics for hosting 
him during the final stages of the writing of this paper. 
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Fig. 1. — Curves for the deep (left) and shallow (right) synthetic surveys showing the 
probability <&/(/#) that a star with magnitude VJ, will be blended with one or more physically 
unassociated stars that have total flux If the primary star has magnitude outside the 
formal magnitude limits of the survey, it must be blended with at least one other star with 
magnitude inside the formal limits (dashed lines). For primary stars with magnitude inside 
the formal limits, there is no such restriction (solid lines). 



Table 1. Properties of the deep and shallow synthetic surveys 





Deep 


Shallow 


FOV longitude (I) 


331.5° 


331.5° 


FOV latitude (6) 


11.0° 


11.0° 


Formal mag. range 


15 < V < 19 


10 < V < 13 


Extended mag. range 


15 < V < 22 


10 < V < 16 


Confusion radius (r c ) 


2" 


20" 


Maximum period (P) 


10 d 


10 d 


Parameters for photometric 






scatter function! : 






a 


+0.018 


+0.029 


b 


-0.357 


-0.477 


c 


-1.161 


-0.726 



tlog 10 a = aV 2 + bV + c (see JO]! 
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Fig. 2. — Fraction of false detections (FDs) in the deep (top left) and shallow (bottom left) 
surveys, and Hot Jupiter detections (HJDs) in the deep (top right) and shallow (bottom right) 
surveys, binned according to the spectral type of the primary. Ten spectral bins are used. 
Results for blended (dotted) and unblended (dashed) configurations are plotted separately, 
as well as combined (solid line). We only plot results for the limiting case that all tertiary 
components are resolved, as no significant qualitative or quantitative change is observed in 
the converse limiting case (see discussion in §4.ip . Results are shown for S/N> 20. 
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Fig. 3. — As in Fig. |2l but for the deep survey only and with detection fractions binned 
according to the effective magnitude, using a bin width of 0.5 mag. Results are shown 
for S/N> 10 and S/N> 30 to illustrate how increasing the SN cutoff biases the detections 
towards brighter magnitudes (see text for discussion). 




Fig. 4. — As in Fig. [3j but for the shallow survey. 
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f B (%) h (%) 

Fig. 5. — Detection fractions for blended configurations, binned according to the percentage 
of the total flux that is due to blended stars, using a bin width of 10%. 



Table 2. Probability p n of n blended stars 





Deep survey 


Shallow survey 


n 


Pn 


Pn 





0.818 


0.528 


1 


0.164 


0.337 


2 


0.017 


0.108 


3 


0.001 


0.023 


4 


0.000 


0.004 
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Fig. 6. — Variation of unblended false detection rates in the deep survey (left) and in the 
shallow survey (right), with respect to transit latitude and the relative size of the eclips- 
ing components, R s /R p . On the grayscale colormap, white corresponds to relatively high 
detection rates and black corresponds to relatively low detection rates. Results are shown 
for S/N> 20 and are effectively identical for both surveys, with most detections caused by 
high-latitude transits. In particular, twins (R p ~ R s ) undergoing high-latitude eclipses are 
the most populous category of false detections. These configurations will produce indistin- 
guishable primary and secondary transits that are shallow enough to be mistaken for a Hot 
Jupiter transit. 
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Fig. 7. — Variation of blended false detection rates in the deep survey (left) and in the 
shallow survey (right), with respect to transit latitude and the blended flux expressed as 
a percentage of the total flux. Grayscale is the same as in Fig. [6] and results are again 
shown for S/N> 20. The most populous category of false detections is those for which the 
blended flux contributes more than 80% of the total flux. A smaller peak is also observed 
for detections in which the blended flux contributes less than 20% of the total flux and the 
transit is high-latitude. 

Table 3. Estimated rates of false detections and Hot Jupiter detections per 10,000 
resolved stars monitored assuming that all tertiary components are resolved* 



Deep survey 





False detections: 




Hot Jupiter 


detections: 






S/N> 


Blended Unblended 


Total 


Blended 


Unblended 


Total 


False/Planet ratio 


10 


2.0 (2.4) 1.8 (1.9) 


3.8 (4.3) 


0.55 (0.69) 


2.7 (3.4) 


3.3 (4.1) 


1.2 (1.1) 


20 


1.4 (1.6) 1.4 (1.6) 


2.8 (3.2) 


0.26 (0.32) 


1.3 (1.6) 


1.6 (1.9) 


1.8 (1.7) 


30 


0.96 (1.1) 1.0 (1.1) 


2.0 (2.2) 


0.14 (0.17) 


0.71 (0.85) 


0.85 (1.0) 


2.3 (2.2) 


Shallow survey 




False detections: 




Hot Jupiter 


detections: 






S/N> 


Blended Unblended 


Total 


Blended 


Unblended 


Total 


False/Planet ratio: 


10 


3.6 (4.4) 0.33 (0.37) 


3.9 (4.8) 


0.30 (0.37) 


0.34 (0.42) 


0.64 (0.79) 


6.1 (6.0) 


20 


3.0 (3.6) 0.34 (0.38) 


3.3 (4.0) 


0.12 (0.14) 


0.14 (0.17) 


0.26 (0.31) 


13 (13) 


30 


2.5 (3.0) 0.32 (0.36) 


2.8 (3.4) 


0.054 (0.063) 0.064 (0.076) 


0.12 (0.14) 


24 (24) 



* Values without parentheses give the predicted detection rates with windowing effects accounted for and the requirement 
that at least three transits are observed, while values in parentheses give the detection rates ignoring any windowing effects. 
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Table 4. Estimated rates of false detections per 10,000 resolved stars monitored assuming 

that all tertiary components are unresolved* 



Deep survey 





Doubles: 




Triples: 








S/N> 


Blended 


Unblended 


Blended 


Unblended 


Total 


False/Planet ratio 


10 


0.34 (0.51) 


0.24 (0.32) 


1.6 (1.9) 


1.5 (1.6) 


3.7 (4.2) 


1.1 (1.1) 


20 


0.23 (0.34) 


0.21 (0.27) 


1.2 (1.3) 


1.3 (1.3) 


2.8 (3.3) 


1.9 (1.7) 


30 


0.15 (0.21) 


0.14 (0.19) 


0.85 (0.93) 


0.96 (1.0) 


2.1 (2.4) 


2.5 (2.3) 


Shallow survey 




Doubles: 




Triples: 








S/N> 


Blended 


Unblended 


Blended 


Unblended 


Total 


False/Planet ratio 


10 


0.66 (0.94) 


0.048 (0.061) 


2.9 (3.3) 


0.25 (0.27) 


3.9 (4.6) 


6.0 (5.8) 


20 


0.53 (0.75) 


0.051 (0.064) 


2.5 (2.8) 


0.26 (0.28) 


3.3 (3.9) 


13 (13) 


30 


0.44 (0.61) 


0.051 (0.065) 


2.1 (2.3) 


0.27 (0.29) 


2.9 (3.3) 


24 (24) 



* Rates presented as in Table [3] Hot Jupiter detection rates are not shown as they do not change from 
Table [3] 



